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Technical Report 

The AM Herculis (hereafter AM Her) objects are short orbital 
per.-'O'^ (P^hours), low mass (M < 2 M^) X-ray binary systems 

7 

composed of a strongly magnetic (B > 10 G) white d^arf and a 
Roche-lobe filling late-type main sequence star. They are the 

strongest known white dwarf X-ray sources having X-ray luminosities of 

3 2 3 3 “I * 

10 - 10 ergs s and bolometric luminosities of 

34 1 

lO"^ ergs s”'*’ . Their X-ray spectra are composed of two 
components. The soft X-ray (E < 2 keV) portions resemble blackbodies 
of temperature ^ 10-50 eV while the hard X-ray (E > 2 keV) 

portions resemble optically thin bremsstrahlung spectra of temperature 
kT^ ^ 10 keV. For a complete review of the observational 
properties of the AM Her objects see Liebert and Stockman (1984) . For 
this grant, we have conducted several theoretical investigations 
concerned with the time variability of AM Her objects, and we have 
endeavored to use HEAO 1-Al observations of these X-ray sources to 
test our theoreticul models. 

Physically, the nature of the AM Her objects is determined by the 
dwarf's large B-field. The B-field ..'synchronizes the systems, prevents 
the formation of accretion disks, and forces field-aligned flows onto 
a fraction f << 1 of the dwarf's surface near its magnetic pole(s). 

The X-ray emission occurs when the funneled plasm^* strikes the dwarf 
and shocks, heating to temperatures of tens of keV. As the 
shock-heated plasma settles onto the dwarf, it cools by 
bremsstrahlung, optically thick thermal cyclotron emission, and 
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Compton cooling. The region in which the plaeraa cools and settles 
onto the dwarf is termed the "emission region." The dominant loss 
process in the emission region is determined by the white dwarf mass 
M*, the "effective" accretion rate (M/f) , and the B-field, 
Bremsstrahlung produces the hard X-rays, while the blackbody limited 
cyclotron emission peaks in the optical or EUV. It is generally 
bel;.eved that the AM Her objects are bremsstrahlung dominated. The 
soft X-ray biackbody component is not produced in the emission region. 
It cones from the dwarf's surfac:e being produced by the reprocessing 
of radiation from the emission region or is simply the intrinsic 
luminosity of the dwarf. For a more complete review of the theoretical 
properties of the AM Her objects see Lamb (1984) . 

Despite extensive research by many workers, the AM Her objects 
are still rather enigmatic. For example, when the observations c2 its 
most well observed member, AM Her, are compared to theory, the 
agreement is rather poor. The largest problem concerns the relative 
strengths of the hard and soft X-ray components. The "standard" 
theory assumes that the blackbody radiation is due to reprocessed 
bremsstrahlung, and so the ratio of the scft X-ray luminosity (L^) 
to the hard X-ray luminosity (Lj^) should be of order 1. In AM 
Her, the ratio is estimated to be anywhere from 5-100. The following 
explanations for this have been advanced: i) the electron scattering 

optical depth through the preshock flow is greater than 1 which 
implies that the bremsstrahlung spectrum is degraded as it traverses 
the preshock flow lowering L^; ii) cyclotron emission is not 
negligible and so lowers the efficiency of bremsstrahlung emission, 
thus lowering iii) electron thermal conduction transports 

energy from the shocked plasma into 'the dwarf lowering and iv) 
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the accreted plasma undergoes steady nuclear burning in the dwarf's 
envelope enhancing L . 

5 

We have investigated the plausibility of items i) , iii) , and iv) 
by performing detailed steady state calculations of the emission 
region structures and of the effects of various preshock flow 
geometries on the emergent spectra. In calculations prior to ours 
(Aizu 1973; Katz 1977; Kylaf is and Lamb 1979) , spherically symmetric 
(SS) flows, equal electron and ion temperatures (a one-temperature 
approach), and negligible electron thermal conduction were assumed. 

We have refined these calculations in two ways: 1) in SS, we allow 

the ion and electron temperatures to differ (a two-temperature 
approach) , include electron thermal conduction, and allow the accreted 
plasma to burn steadily (see Imamura et al . 1984a; and Weast et al . 
1984) ; and 2) in a one-temperature approach v,’ithout electron tnermal 
conduction, we relax the assumption of SS by considering dipolar field 
aligned flow (see Imamura and Durisen 1983) . Our work on item 1) with 
D. Q. Lamb and G. J. Weast was initiated in 1978. Although net 
directly funded by the grant, it has been continued in parallel with 
our funded research. Because of its usefulness to the grant-related 
investigations, we summarize it in the next two paragraphs. Our work 
on item 2) was undertaken specifically for the purposes of this grant 
and is described in subsequent paragraphs. 

From an observational viewpoint, two-temperature and electron 
thermal conduction effects do not significantly alter the properties 
of accreting dwarfs over most of parameter space if the accreted 
plasma does not undergo steady nuclear burning. Only for M* > 1.3 
is a two-temperature approach required to calculate accurate 
luminosities and X-ray temperatures. In these cases, Lj^ is 



Page 5 


increased by up to a factor of 1.5 and kT^ is decreased by up to a 

factor of 2. Concerning the emission region structures, however, 

two-temperature and electron thermaJ conduction effects are 

significant over a large portion of parameter space. For M* > 1.0 

a two-temperature approach with electron thermal conduction is 

required to calculate accurate emission region structures. We note 

that two-temperature and electron thermal conduction effects are 

important whenever Compton cooling is important. Further, we note 

that electron thermal conduction only transports energy around the 

emission region and does not act as a global loss mechanism. This 

means that electron thermal conduction is not the reason why 

L/L. >> 1. 
s h 

If the accreted plasma is allowed to undergo steady nuclear 
burning, the blackbody luminosity of the dwarf is greatly increased. 
Typically Lj^j^ increases by the factor 

(QXM)/(GM*M/R*) = 17 (R*/5xl0®cm) 


where the denominator is the energy gained by the plasma falling into 
the dwarf's gravitational potential well and the numerator is the 
energy released by nuclear burning of the accreted plasma. The 
individual terms are G, the gravitational constant, Q, the energy 
released per gram in the proton-proton process, R*, the white 
dwarf radius, and X, the hydrogen mass fraction of the plasma. The 


enhanced 


increases 


the -adiation field in the emission region 


which enhances the effects of Compton cooling and therefore. 


two-temperature effects. When the accretec. pi*sma undergoes steady 


nuclear burning, a two-temperature approach is required for an 
accurate treatment at a] 1 masses of interest. The major effects of 
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nuclear burning are to depress and kT^ and to enhance 
L , The depressed kT makes nuclear burning models 
unattractive solutions to the AM Her dilemma because AM Her has 
kT^ = 31 keV, with an uncertainty of about 5 keV, while the 
hardest a nuclear burning dwarf can appear is about 26 keV (at M^ 

= 1.4 M^) . Thus, AM Her is only marginally consistent with ':hese 
calculations. 

F: eld-aligned flows tend to produce X-ray spectra which are quite 
similar to ones produced in SS calculations. However, the "funneling" 
does le, d to: i) radiation escaping to infinity with a "fan beam" 
emission pattern; and ii) smaller electron scattering optical depths 
to infinity ( T-^) when compared to SS ones of the same effective 
accretion rate which leads to less degradation of the bremsstrahlung 
spectra. Concerning item i) , we have used the "fan beam" emission 
patterns to calculate light curves magnetic white dwarfs in binary 
systems are likely to produce for a wide range of orbital inclinations 
and magnetic field orientations. The light curves are rather 
sensitive functions of these parameters and can be used to place 
constraints upon the geometry of AM Her's emission region. We are 
currently in the process of comparing our light curves to the observed 
light curves of AM Her. We had proposed to use HEAO 1-Al data to 
generate additional light curves but were disappointed to find, after 
considerable effort, that the data were too noisy. Our light curves 
did not improve on others in the literature. 

Using our theoretical models, we have also considered the 
question of the large ratio (Imamura 1984) . The 

calculations suggest that degradation of the bremsstrahlung can't 
explain the large ratio in AM Her. To get the observed 
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ratios, requires that the degradation be very severe. The large 
degradation requires a large which has two problems: 1) it 

leads to kT^ << 30 keV; and 2) it implies that (or 

equivalently f << 1) , which leads to >> 50 eV, Here Mg 

(«41TR^c/^g) is the Eddington accretion rate. Both of these 
results do not agree with the observations. 

The previous discussion has assumed that a steady state model of 
the AM Her objects is reasonable. However, recent work, both 
theoretical and observational, suggests that AM Her objects may be 
variable. Observations of the AM Her objects AN UMa and El 405-451 
have shown that they show coherent periodicities on time scales 
ranging from 1 to 2 seconds in their optical emission (Middleditch 
1982) . Middleditch suggests that these variations are due to the 
oscillatory instability of the emission region structures found by 
Langer, Chanmugam- and Shaviv (1981,1982) and Langer (1984) during 
time- dependent calculations of accretion flows onto magnetic white 
dwarfs. Langer et al. found that one-temperature shocks with power 
law cooling functions proportional to p^T** are unstable to 
oscillatory motions whenever oi< 0.6. Because AM Her shocks are 
probably bremsstrahlung dominated (i.e. = 0.5), Langer et al. 

suggested that they are likely candidates to show this instability. 
They predict that the oscillations will occur with time scales on the 
order of three times the postshock plasma cooling time scale 
(tcooi^ , which ranges from tenths of a second to several seconds, 
consistent with the observations of Middleditch. However, because 
Middleditch observed the periodicities in the optical it is not clear 
that the variations reflect variations in the emission region 
structure. The optical emission is not produced in the emission 
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region. It probably arises from a region higher up in the preshock 
funnel. 

To further the understanding of the nature of the periodicities 

seen by Kiddleditch, we have also studied the stability properties of 

the emission regions. Our calculations differ from Langer et al.'s in 

that we consider SS flowsr and we use a two- temperature approach with 

electron thermal conduction and Compton cooling in some of our 

calculations (see Imamura et al. 1984b, c) . With regard to the 

different geometries, the differences are small because the 

thicknesses of the emission regions are always much less than R*. 

Our results qualitatively agree with Langer et al., but differ in 

detail. We find that: i) several oscillatory modes are possible, 

which we call in order of increasing oscillation frequency the 

Fundamental (F) , the First Overtone (10) , the Second Overtone (20) , 

and so on; ii) the F, 10, and 2'' are unstable for oi < 0.4, 0.8, and 

0.8; iii) the F, 10, ana 20 have oscillation periods of 

5/3t^^Q^, and iv) for bremsstrahlung dominated shocks 

(i.e., cK- 0.5), there will be periodicities on the time scales 

and fractional luminosity variations of up 

to 50%; and v) for accretion onto a 1.0 M white dwarf, the 

o 

inclusion of two-temperature effects, electron thermal conduction, and 
Compton cooling stabilizes all oscillatory modes. The results of our 
nonlinear calculations are in excellent agreement with the linear 
analysis of Chevalier and Imamura (1982) . Langer et al. appear to see 
the F but no overtones. We do not have an explanation for this 
discrepancy. 

Observing this instability would be quite valuable because 
^cool' t)eing approximately given by 
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tcooi' « 0.012(Mgf/M) seconds, 

is a function of the effective accretion rate and the depth of the 
dwarf's gravitational potential well. In the estimate for 
I is the specific internal energy, is the bremsstrahlung 

volume emissivity, and we have used the Rankinc-Hugoniot strong shock 
jump conditions to define typical postshock plasma conditions. Thus, 
if M/f is known, a direct determination of the dwarf's gravitational 
potential can be made. Further, if one assumes a white dwarf 

4 

mass-radius relationship (see Chandrasekhar 1957} , the mass and radius 
of the dwarf can be separately determined. 

We are currently in the process of analyzing pointed data of AM 
Her taken by the HEAO 1-Al experiment. The A1 experiment was 
sensitive to photons in the energy range 0.5-20 keV and had a temporal 
resolution of 5 ms, so in principle, variations in the emission region 
structure would be expected to show up in the data. However, because 
AM Her shows substantial "flickering" on time scales of tens of 
seconds and the signal to noise ratio (S/N) of the time-averaged data 
is /vl, the detection of low amplitude, short period coherent 
fluctuations is difficult. We are currently in the process of 
searching for periodicities in the A1 data using Fourier technigiues. 

To summarize, we have made several useful and important 
contributions to our theoretical understanding of X-ray variability in 
AM Her objects. We have generated funneled accretion models that make 
detailed predictions of X-ray light curves (Imamura and Durisen 1983) . 
When applied to real systems, they constrain various geometrical and 
physical parameters (Imamura 1984) , The same funneled accretion 
models have provided useful tests of various explanations for the 


large observed (Iinamura 1984) . We have also 

investigated the oscillatory stability of the postshock regions of 
white dwarf accretion flows (Imamura et al. 1984b, c). Unfortunately, 
the BEAO l-Al scan data did not yield quality light curves for 
comparison with our funneled accretion models, but we are now in the 
process of Fourier analyzing pointed data in the hopes of detecting 
evidence for shock oscillations. 
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We also plan to have a publication on the Fourier analysis of 
HEAO 1-Al data. 



